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Abstract 

Localization methods reduce the path integrals in Af > 2 supersymmetric Chern-Simons 
gauge theories on S 3 to multi-matrix integrals. A recent evaluation of such a two-matrix 
integral for the M = 6 superconformal U (N) x U (N) AB JM theory produced detailed 
agreement with the AdS/CFT correspondence, explaining, in particular, the N 3 ^ 2 scaling of 
the free energy. We study a class of p-matrix integrals describing M = 3 superconformal 
U (N) p Chern-Simons gauge theories. We present a simple method that allows us to evaluate 
the eigenvalue densities and the free energies in the large N limit keeping the Chern-Simons 
levels ki fixed. The dual M-theory backgrounds are AdS^ x Y, where Y are seven-dimensional 
tri-Sasaki Einstein spaces specified by the ki. The gravitational free energy scales inversely 
with the square root of the volume of Y. We find a general formula for the p-matrix free 
energies that agrees with the available results for volumes of the tri-Sasaki Einstein spaces 
Y, thus providing a thorough test of the corresponding AdS 4/ CFT 3 dualities. This formula 
is consistent with the Seiberg duality conjectured for Chern-Simons gauge theories. 
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1 Introduction and Summary 



The AdS/CFT correspondence [lj|3] provides many predictions about the dynamics of strongly 
interacting field theories in various numbers of dimensions. For the case of three dimensions, 
ref. [2] predicted that the number of low-energy degrees of freedom on N coincident M2- 
branes scales as N 3 ^ 2 for large N. Remarkably, this surprising prediction was recently con- 
firmed [5] in the context of the ABJM construction [6] of the U(N)k x U(N)-k Chern-Simons 
gauge theory on coincident M2-branes. The paper [5] was in turn based on [7] where the 
methods of localization [8] were shown to reduce the path integral of the Euclidean ABJM 
theory on S 3 to a matrix integral. This matrix model was solved in the large iV limit with 
N/k kept fixed [5|[9], leading to precise tests of the AdS/CFT correspondence for Wilson 
loops and for the free energy (by which we mean minus the logarithm of the Euclidean par- 
tition function). The exact solution of this matrix model is related by analytic continuation 
to a solution [lO] of another matrix model describing topological Chern-Simons theory on 
S 3 /^; in particular, the formula for the resolvent has the same structure in the two cases. 
A generalization of the matrix model to the case where the Chern-Simons levels do not add 



up to zero was considered in 11 



The aim of our paper is to build on the major progress recently achieved in [5j|7||9] in 
several ways. In section 2 we revisit the matrix integral for the ABJM theory on S 3 and 
uncover the details of the eigenvalue distribution. The matrix eigenvalues are located along 



the branch cuts of the resolvent used in [5] and derived in 10 for the 5 ,3 /Z 2 model. While 
the endpoints of the cuts can be read off directly from the resolvent, the cuts themselves are 



not simply parallel to the real axis, in contrast with the matrix model of 10 . In order to 
gain intuition for the location of the eigenvalues, we develop a numerical method for finite 
N and k. This method allows us to access values of N and k that are large enough for the 
result to be a good approximation to the limit studied in [5j|9]. Furthermore, we focus on 
the limit where N is sent to infinity at fixed k where the ABJM model is expected to be 
dual to the AdS^ x S 7 /Z& background of M-theory. In this strong coupling limit, which is 
not of the 't Hooft type, we find analytically that the structure of the solution simplifies 
considerably. An ansatz where the real parts of the eigenvalues scale with allows us to 
calculate the free energy analytically. Unlike in [H], our method does not rely on resolvents 
or mirror symmetry. We confirm that the free energy scales as N 3 / 2 with the coefficient 
found in [H]. 

In section 3 we develop our analytic approach further and apply it to the large N limit 
of matrix models describing quiver Chern-Simons gauge theories on S 3 . We study explicitly 
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a class of M = 3 superconformal U(N) P gauge theories with bifundamental matter, quartic 
superpotentials, and Chern-Simons levels k\, k 2 , • • • , k p that sum to zero. These models were 



introduced in 12 where their type IIB brane constructions were presented. The type IIB 



brane constructions involve iV D3-branes that are wrapped around a circle and intersect the 
(1, qi), (1, q 2 ), . . . , (1, q p ) 5-branes located sequentially along the circle. The dual AdS^ x Y 
M-theory backgrounds for these models, which involve certain seven-dimensional tri-Sasaki 



Einstein spaces Y, were conjectured in 13 . The tri-Sasaki Einstein spaces Y are defined 



to be bases of hyper-Kahler cones 14 16 , and we take the Einstein metric on them to be 
normalized so that R mn = 6g mn . The p-matrix models for the gauge theories dual to AdS^ x Y 
may be read off from [7] . In the large N limit we calculate the eigenvalue densities for these 
matrix models and show that they are piecewise linear. This remarkably simple conclusion 
allows us to evaluate the coefficient of the N 3 ^ 2 scaling of the free energy as a function 
of the levels ki and compare it with the calculation on the gravity side of the AdS/CFT 
correspondence [5 17 . For an arbitrary compact space Y we find that the gravitational free 
energy is 



TV 3 / 2 



2tt 6 



27Vol(F) 



For p = 3 the tri-Sasaki Einstein spaces Y are the Eschenburg spaces 18 whose volumes 



were determined explicitly in 19 . Our 3-matrix model free energy is in perfect agreement 
with this volume formula. 

Furthermore, we carry out calculations of the p-matrix model free energy and use them 
to conjecture an explicit general formula for the volumes via the AdS/CFT correspondence. 
For a general p-node quiver with CS levels k a = q a+x — q a , with 1 < a < p and q p+ \ = qi, we 
conjecture in section 4 that 



Vo\(Y) _ E(v,g)er H(a,b)eE Ife ~ g&l 
Vo\{S 7 ) " 



6=1 \Qa - qb\ 



(2) 



where the sum in the numerator is over the set T of all trees (acyclic connected graphs) 
with p nodes. Such a tree (V, E) consists of the vertices V = {1,2, ... ,p} and \E\ = p — 1 
edges. The volumes of the corresponding tri-Sasaki Einstein spaces Y had previously been 
studied by Yee, who expressed them through an integral formula (eq. (3.49) of [20]). In 
the cases we have checked, our formula ^ is consistent with that of [20]. Equation ^ is 
invariant under permutations of the q a , supporting the conjectured Seiberg duality for Chern- 
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Simons theories with at least M = 2 supersymmetry 21-23 , which may be motivated by 
interchanging different types of 5-branes in the type IIB brane constructions of these models. 



Recent work [5j|7\24 and the present paper hint at a special role that may be played 
by the Euclidean path integral of a three-dimensional conformal field theory on S 3 . This 
quantity may be analogous to the conformal anomaly coefficients in even dimensional CFTs. 
Recall that the anomaly coefficients are very useful measures of the number of degrees 
of freedom. For even dimensional theories with weakly curved dual backgrounds, these 



coefficients can be calculated using dual gravity in AdS space 25 leading to precise tests 
of the AdS/CFT correspondence. Such a definition of the number of degrees of freedom is 
not available for three-dimensional CFTs. As mentioned already, the path integral on S 3 
can be reduced to matrix integrals using supersymmetric localization methods [71. Earlier 



work on gravity in Euclidean AdS^ 17 has pointed to the usefulness of the corresponding 
quantity: After adding certain surface counter-terms, the action becomes finite and appears 
to be unambiguous. The successful matching of this finite gravitational action with the path 
integral on S 3 in [5] for the M = 6 ABJM theory and in the present paper for a class of 
J\f = 3 superconformal theories provides evidence for the usefulness of this quantity as a 
measure of the number of degrees of freedom. 

One may hope that the free energy on S 3 is also a useful quantity for nonsupersymmetric 
fixed points]^] For example, one could aim to match the large N free energy on S 3 for the 
nonsupersymmetric example of AdS/CFT correspondence conjectured for the O(N) sigma 
model in three dimensions |26|. 



2 ABJM Matrix Model 
2.1 Matrix Model Setup 

As shown in [71, the partition function for ABJM theory on S 3 localizes on configurations 
where the auxiliary scalars a and a in the two M = 2 vector multiplets are constant N x N 
Hermitian matrices. Denoting the eigenvalues of a and a by Aj and Aj, with 1 < i < N, one 

1 Of course, another nonsupersymmetric measure of the number of degrees of freedom, which is very useful 
physically, is the thermal free energy. 
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can write the partition function as 



n d\j d\j \ (2sinh^) ( 2 sinh^ (ik^ iV2 y2 



(3) 



where k is the Chern-Simons level, and the precise normalization was chosen as in [5j. The 
integration contour should be taken to be the real axis in each integration variable. When 
the number N of eigenvalues is large, the integral in eq. ^ can be approximated in the 
saddle-point limit by Z = e~ F , where the "free energy" F is an extremum of 



F(A il A i ) = -iA^( A ;-Aj)-X;iog 



47T 

3 KJ 



2sinh^^| [2sinh^^ 



2 

Aj — Xj 



(4) 



+ 2 log 2 cosh - J ' + 2 log AM + 2iV log(27r) 



i,3 



with respect to Aj and Aj. The goal of this section is to compute the leading contribution to 
F in such a large iV expansion while holding k fixed. 

Varying (Hi) with respect to Aj and Xj we obtain the saddle-point equations: 



dF ik \ - Xj-Xi , v-^, Xj-Xi 



Xi — coth — + tanh 



dXi 2ti 1 ^ 2 ^ 2 

dF ik~ ^ Xj - \ , v-^ , , Aj - A 



(5) 



<9Aj 2tt 



Aj - coth j i + tanh j i = . 



Similar saddle-point equations appear in the context of matrix models derived from topolog- 



ical string theory 10 27 , and can be solved using powerful techniques based on holomorphy 
arguments and special geometry. Such methods were used in |5j to solve eqs. ^ in the limit 
where N is taken to infinity while holding N/ k fixed. Our goal in this paper is more modest 
than solving the matrix model for any value of the 't Hooft parameter N/k. We will work 
at fixed k and take the large N/k limit. As we will see shortly, in such a limit we can find 
the eigenvalue distribution using a more elementary approach. 
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2.2 A Numerical Solution 



To gain intuition, one can start by solving the saddle-point equations ^ numerically for 
any values of iV and k. One of the simplest ways to do so is to view equations ^ as 
describing the equilibrium configuration of 2N point particles whose 2-d coordinates are 
given by the complex numbers Aj and Aj and that interact with the forces given by eq. ([5]). 
This equilibrium configuration can be found by introducing a time dimension and writing 
down equations of motion for \j(t) and \j(t) whose solution approaches the equilibrium 
configuration (|5| at late times in the same way as the solution to the heat equation approaches 
a solution to the Laplace equation at late times. The equations of motion for the eigenvalues 
are 

d\i _ OF d\ _dF 

where t\ and t~ x are complex numbers that need to be chosen in such a way that the saddle 
point we wish to find is an attractive fixed point as t — > oo. 

N=25 and 100, k=l 
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Figure 1: Numerical saddle points for the ABJM matrix model. The eigenvalues for N = 25 
are plotted in black and those for N = 100 are plotted in orange. The plot has been obtained 
with r\ = tt 



1. As mentioned in the text, the real parts of the eigenvalues grow with y N 



In figure [T] we show typical eigenvalue distributions that can be found using the method 
we just explained. There are several features of the saddle-point configurations that are 
worth emphasizing: 

• The eigenvalues Aj and A, that solve ^ are not real. 

That the eigenvalue distributions do not lie on the real axis might be a bit puzzling 
given that Aj and Aj are supposed to be eigenvalues of Hermitian matrices. However, 
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it is well known that in general, when using the saddle-point approximation, the main 
contribution to an oscillatory integral may come from saddles that are not on the 
original integration contour but through which the integration contour can be made 
to pass. We will assume that the integration contour that should be chosen in writing 
down the integral in eq. ^ can be deformed so that saddle points like those in figure [j] 
are the only important ones. 

• The eigenvalue distributions are invariant under Aj — > — Aj and Aj — > — Aj. 

Indeed, the saddle-point equations ^ are invariant under these transformations, so it 
is reasonable to expect that there should be solutions that are also invariant. 

• In the equilibrium configuration the two types of eigenvalues are complex conjugates 
of each other: Aj = Xj. 

Indeed, it is not hard to see that upon setting Aj = Aj the two equations in (|5| become 
equivalent, so it is consistent to look for solutions that have this property. 

• As one increases N at fixed k, the imaginary part of the eigenvalues stays bounded 
between — 7r/2 and tt/2, while the real part grows with N. We will show shortly that, 
for the saddle points we find, the real part grows as N 1 ^ 2 as iV — > oo. 



2.3 Large N Analytical Approximation 

Let us now find analytically the solution to the saddle-point equations ^ in the large N 
limit. As explained above, we can assume Aj = Aj and writ^] 

Aj = N a Xj + iyj , Aj = N a Xj — iyj , (7) 

where we introduced a factor of N a multiplying the real part because we want Xj and y 3 - to 
be of order 0(N°) and become dense in the large N limit. The constant a is so far arbitrary 
but will be determined later. 

In passing to the continuum limit, we define the functions x, y : [0, 1] — > R so that 

x 3 = x(j/N) , y 3 = y{j/N) . (8) 
2 After completing this work, we became aware that ref. |28| employs a similar ansatz. 
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Let us assume we order the eigenvalues in such a way that x is a strictly increasing function 
on [0, 1] . Introducing the density of the real part of the eigenvalues 



P{x) = £ , (9) 
dx 



one can approximate (El) as (see Appendix [A|) 

F = -N 1+a [ dx xp{x)y(x) + N 2 ~ a I dx p(x) 2 f(2y(x)) + ■ ■ ■ , (10) 



7T 



where the function / is 

/(t) = 7r 2 -(arge t *) 2 . (11) 

In other words, / is a periodic function with period ir given by 

/(*) = vr 2 - t 2 when - tt < t < tt . (12) 

It may be a little puzzling that while the discrete expression for the free energy in eq. Q 
is nonlocal, in the sense that there are long-range forces between the eigenvalues, its large A" 



limit (10) is manifestly local. One can understand this major simplification from examining, 
for instance, the first saddle-point equation in The force felt by Aj due to interactions 
with far-away eigenvalues Aj and Aj is 

— coth ^ + tanh ^ « — sgn(Re Aj — Re Aj) + sgn(Re Aj — Re Aj) , (13) 

the corrections to this formula being exponentially suppressed in Re Aj — Re A, and Re Aj — 
ReAj. In other words, the nonlocal part of the interaction force between eigenvalues is 



given just by the right-hand side of eq. (13). The nonlocal part of the force vanishes when 
ReAj = ReAj, so in assuming that the two eigenvalue distributions are complex conjugates 
of each other, we effectively arranged for an exact cancellation of nonlocal effects. All that is 



left are short-range forces, which in the large limit are described by the local action (10). 



One can view F as a functional of p(x) and y(x) and look for its saddle points in the set 



C = |(p, y) : J dx p(x) = l',p(x) > pointwisej- 



(14) 
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These constraints mean that p is a normalized density. Motivated by the numerical analysis 
we performed, we assume that p and y describe a connected distribution of eigenvalues 
contained in a bounded region of the complex plane. 



Let us assume a saddle point for F exists. As N — > oo, we need the two terms in (10) to 



be of the same order in N in order to have nontrivial solutions, so from now on we will set 



a 



1 
2 



(15) 



The real part of the eigenvalues therefore grows as iV 1//2 , and to leading order, the free energy 
behaves as N 3 ^ 2 at large N. In writing (10) we omitted the last two terms from eq. Q. 
They do not depend on p or y and hence do not affect the saddle-point equations. They are 
also lower order in N given the choice of a. 



To find a saddle point for F, one can add a Lagrange multiplier p to (10) and extremize 



p = N 3/2 



— [ dxxp(x)y(x) + [ dx p(x) 2 f(2y(x)) — — 

7T ./ ./ 2lX 



dx p(x) 



(16) 



instead of (10). As long as p{x) > 0, the saddle-point eigenvalue distribution satisfies the 
equations 



4:irp(x)f(2y(x)) = p — 2kxy(x) 
2np(x)f'(2y(x)) = -kx . 



(17) 



Plugging (12) into (17) one obtains 

p{x) 



_p_ 
4tt 3 ' 



y(x) 



TT 2 kx 

2p 



[IS) 



as long as — n/2 < y(x) < tt/2. If p is supported on [— for some > that we will 
determine shortly, we can calculate p from the normalization of the density p(x): 



dx p(x) = 1 



2tt 3 



(19) 



Plugging this formula back into (fTol), we obtain the free energy in terms of x* 



iV 3 / 2 (127T 4 + k-X 

24ir 2 x* 



+ o(iV 3 / 2 ) 



(20) 



S 



This expression is extremized when 

vM = \. (21) 

Luckily, the answer y(x*) = 7r/2 is consistent with our assumption that — 7r/2 < y(x) < 



7T 



7r/2 without which eq. (18) would not be correct. In Appendix |B| we check that assuming 



y(x*) > 7r/2 implies a contradiction. The extremum of F obtained from eqs. (20) and (21) 
is 



F = ^k^N 3 / 2 + o{N 3 ' 2 ) . 



(22) 



This result agrees with the free energy found in [5] using the regularized Euclidean action in 



AdS, 17 
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Figure 2: Comparison between analytical prediction and numerical results for the density of 
eigenvalues p defined in eq. The dotted black line represents the analytical calculation, 
and the numerical result is shown in orange dots. 

In the large iV limit the eigenvalues therefore condense on two line segments, and on 
these two line segments they have uniform density. In figure [2] we compare the analytical 
result for the density with the numerical one. 

We would like to compare the location of our eigenvalue distributions with the results 
of |5|. Noting a similarity between the ABJM matrix model and the S 3 /7j2 matrix model 



solved in 10 , Drukker, Marino, and Putrov [5] wrote down a resolvent for the ABJM model. 
This resolvent has cuts in the A plane corresponding to the locations of the eigenvalues. In 
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particular, it has a cut where the Aj eigenvalues are located and a second cut where the Aj 
eigenvalues are located but shifted by ni. More specifically, the resolvent has the form 



u(X) = 2 log ^- |y( e A + &)( e A + i/ 6 ) - y/(e*-a)(e*-l/a) 
where a + l/a + 6+ l/6 = 4 and at strong coupling, 



(23) 



1,1 / 2N 1 



(24) 



The ellipses denote terms exponentially suppressed in N/k relative to the leading term. 
Solving the equations for a and b, we find that the branch points in the A plane are at 



± log a = 7i 



± log b 



-7J-1 



27V 


1 


~k ~ 


12 


/2N 
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~k ~ 


' 12 



ITT 

in 
' ~2 ' 



(25) 
(26) 



These expressions are in agreement with (21) in the large N limit. 



Let us also try to compare our findings with the exact results found for the supersymmet- 
ric Wilson loops in ABJM theory [5j[9]. The expectation values of 1/6 and 1/2 supersym- 
metric Wilson loops are proportional, respectively, to the expectation values of J2f=i &Xi an d 
Y^iLi eXi + eXl i n the matrix model 5,7,9,29 . In our approach, these quantities become 



<< 6 > 



2niN 



e x{x) p(x)dx . 



(W^) = ^ r (e*W + e*<»> ') p(x)dx . 



(27) 
(28) 



If we evaluate (27) and (28) using the saddle point we have found, we get 



(W# 8 > 




(29) 

(30) 



The exponents in these formulae agree with the results in [5j|9 29 . 

We should keep in mind that the ABJM model has a type IIA string interpretation only 
in the limit where N/k ^> 1, N 1 ^ 2 /A; 5 / 2 <C 1. These conditions apply only in the limit 
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Figure 3: Necklace quiver diagrams for U(N) P Chern-Simons gauge theories. 

where both N and k are taken to infinity. Our approximations are only applicable in the 
M-theory limit where N is taken to infinity at fixed k. Thus our Wilson loops have a dual 
interpretation as wrapped M2-branes in M-theory rather than as strings in type IIA string 
theory. 

3 Necklace Quiver Gauge Theories 

In this section we consider a class of quiver Chern-Simons U{N)k 1 x U(N)k 2 x • ■ ■ x U (N)k p 
gauge theories whose quiver diagrams look like necklaces (see figure [3J [l2j[l3] . In the M — 2 
superspace formulation the theory is coupled to bifundamental chiral superfields A a and B a , 
a = 1,2, ... ,p, whose interactions are governed by the quartic superpotential 



W 



2tt 



J2-rr(Ba-iA a ^-A a B a ) 

t\i. 



a=l 



(31) 



The 4-d "parent theories" for these Chern-Simons models, i.e. the 4-d gauge theories with 
the same type of quiver diagrams and superpotentials, arise on iV D3-branes placed on 
generalized conifolds 30,31 . For Y7 a =i^a = and N = 1 the moduli space of each 3-d 



model was calculated in 12, 13 and shown to be given by a certain hyper-Kahler cone in 



four complex dimensions 14-16 . Therefore, it was conjectured that such a model describes 



the low-energy dynamics of N coincident M2-branes placed at the tip of this cone. The 
quiver Chern-Simons gauge theories were therefore conjectured to be dual to the AdS^ x Y 
backgrounds of M-theory where Y are the bases of the cones. These gauge theories are 
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natural generalizations of the ABJM model that is found for p = 2. 

The gauge theories in figure [3] arise from the following type IIB brane construction [12] . 
We consider N D3-branes filling directions 0123 where x 3 is a circle, and then add a sequence 
of (1, q) 5-branes (the bound states of an NS5-brane and q a D5-branes) filling directions 012 
and located at specific points along the x 3 circle. If we express the CS levels as 

K = q a +i ~ Qa , (32) 

where q a are integers, then as we move around the circle we find the ordering (l,q a ) of 
the 5-branes. This picture is helpful in considering certain transformations in the gauge 



theory that are analogous to the Seiberg duality in four-dimensional gauge theory 21-23 
As in that case, these transformations are related to interchange of adjacent branes and thus 
correspond to interchange of q a and q a +i- (There is also a shift in the rank of one of the 
gauge groups that may be neglected in the large N limit.) Our explicit answer for the free 
energy will have this symmetry. 

For a general set of Chern- Simons levels such &p > 2 gauge theory has M = 3 superconfor- 
mal invariance, but in the special case where p is even and the CS levels are (yk ? k ? k ; k ; . . . ) 



the supersymmetry is enhanced to M = 4 12,32,33 . Then Y = S ,7 /(Z p / 2 x Zfc p / 2 ) and 



Vol(y) = 47r 4 / (3kp 2 ) . For more general ki the eight-dimensional cone is not an orbifold, 
which complicates the calculation of its volume. Nevertheless, these volumes were computed 
in 19,20 , and we can compare the result on the gravity side with our calculation of the free 



energy. 



3.1 Multi-Matrix Models 

As explained in [7] , the partition function for the necklace quivers in figure [3] localizes on 
configurations where the scalars a a in the M = 2 vector multiplets are constant Hermitian 
matrices. Denoting by A Q) i, 1 < i < N, the eigenvalues of a a , the partition function takes 
the form of the matrix integral 



n " |nf" :v ' ; > ex P 



Vif 2n :i\ n M 2cosh^ 



,,,,,, 



i 

An 



(33) 



/ 



The normalization of the partition function was chosen so that it agrees with the ABJM 
result from eq. ^ in the case p = 2. As in the ABJM case, the integration contour should 
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be taken to be the real axis in each integration variable. The saddle-point equations following 



from (33) are 

— X a>i - 2 2^ coth — '- + 2^ tanh + ^ tanh = . (34) 

j^i 3 3 



These equations can be solved numerically using the method described in section 2.2| By 



replacing the right-hand side of these equations by r a d\ a j/dt, we obtain a system of first 



order differential equations whose solution converges at late times t to a solution of eq. (34) 
provided that the constants r a are chosen appropriately. We will now show how to obtain 
an approximate analytical solution valid in the limit where N is taken to be large and k is 
held fixed. 

Based on our intuition from the ABJM model, let us assume that in this case too the 
real part of the eigenvalues behaves as N 1 ^ 2 at large N while the imaginary part is of order 
one. So if one writes 

Xaj = N 1/2 x aJ + iy ad , (35) 

then the quantities x a j and y a j become dense in the large N limit. Under this assumption, 
we will be able to solve the saddle-point equations to leading order in N in a self-consistent 
way. We can pass to the continuum limit by considering the normalized densities p a (x) of 
the x a> j together with the continuous functions y a (x) that describe the imaginary parts of the 
eigenvalues as functions of x. Let us first make a rough approximation to the saddle-point 



equations (34). When N is large, we have 



coth ^ a ' 3 ^ a ' 1 w sgn (x a j - x aji ) , tanh ^ ^ a±1 ' 1 ~ sgn ( Xa j _ Xa±l i ) . (36) 

To leading order in N, the saddle-point equations then become 

J dx' [2p a (x') - p a +i{x') - Pa-iW)) sgn(x' - x) = . (37) 

Differentiating with respect to x, we immediately conclude that all p a must be equal to one 
another to leading order in N, so we can write p a (x) = p(x) for some density function p(x) 
that is normalized as 



/ 



dxp{x) = l. (38) 
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With the simplifying assumption that the densities p a are equal, one can go back to the 



integral (33) and calculate the free energy functional F[p,y a ] to leading order in N (see 



Appendix |Aj): 

AT3/2 



N 6 ' 2 f vf^ 
F \P,ya\ = ~2 — / dxxp(x)2_^k a y a (x) 

71 J 0=1 

at 3 / 2 r p 
+ — / dx P (x) 2 Y,f(y a+ i(x)-ya(x)) + o(N 3 / 2 ) 

A J 0=1 



where / is the same function that was defined in (11). We wish to evaluate the integral 



(33) in the saddle-point approximation where it equals Z = e , the free energy F being 
an appropriate critical point of F[p,y a }. Let us assume that the eigenvalue distribution 
corresponding to this saddle point is connected, symmetric about x = y = 0, and bounded. 



In looking for the eigenvalue distribution that extremizes (39) to order 0(N 3 ^ 2 ), an 
important observation is that, in fact, one cannot find this distribution, because to this 
order in N F[p,y a ] has a flat direction given by y a {x) —> y a (x) + 8y(x) for any function 
5y(x). The second term in eq. ( [39] ) is clearly invariant under this transformation, and the 
first term is also invariant because J^o=i ^ a = 0- The existence of this flat direction is not 
a problem at all if one just wants to compute the free energy F to leading order in N. If 
one's goal is instead to find the eigenvalue distributions for the saddle point, subleading 



corrections to (39) that presumably lift this flat direction must be taken into account. In 
this paper we will content ourselves with calculating the free energy to order 0(N 3 ^ 2 ), and 
will leave a careful analysis of how the flat direction gets lifted for future work. 



Before we examine the extremization of the free energy functional (39) in more detail 



let us make a few comments and present a result that follows already from the discussion 



above. Suppose we manage to find a saddle point of F by extremizing (39) for a quiver 
Chern-Simons gauge theory that in the large N limit and at strong 't Hooft coupling is dual 
to an AdS^ x Y M-theory background. Let us assume that this saddle point gives the most 



important contribution to the partition function. What can we learn? From (39) one may 
infer that the free energy grows as iV 3 / 2 at large N as expected from supergravity, so our 
computation provides a gauge theory explanation of this N 3 ^ 2 behavior. Moreover, one can 
compare the free energy we obtain with the exact M-theory result 

F = N 3/2 \ — — — (40) 
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that can be derived as a straightforward generalization of the ABJM computation in p\. Via 
this formula we will compare successfully our matrix model results with the expressions for 
the volumes of tri-Sasaki Einstein space available in the literature 19 , 20 . 



3.2 A Class of Orbifold Chern-Simons Theories 

The vacuum moduli space of the nonchiral quivers with alternating CS levels (k, —k, k, —k, . . .) 
and N = 1 is the orbifold C 4 / (Z p/2 x Z kp/2 ) |12]. There is an induced orbifold action on 
the unit 7-sphere in C 4 , and thus the internal space Y is S 7 / (Z p / 2 x Zkp/2) ■ Consequently, 
we expect 

where in the second equality we used the round 7-sphere volume Vo\(S 7 ) = 7r 4 /3. 

This formula can be reproduced very easily from the matrix model computation. The 



saddle-point equations (34) are solved by setting \ 2a ,i = A, and A2 +i,i = X{, \ an d A, 
being the eigenvalues for the saddle point of the ABJM matrix model discussed in detail in 
section [2] The free energy of the p-node quiver with CS levels (k, —k, k, —k, . . .) is therefore 
p/2 times the free energy in the ABJM model, and thus 

F = |F ABJM = ^V /2 iV 3/2 + o{N^) . (42) 
2 



Using eq. (40), one immediately reproduces the volume of the S 7 orbifold in eq. (41). 



3.3 Warm-up: A Four-Node Quiver 

Another case we can easily solve using the approximation scheme developed above is that 
of the four-node quiver with CS levels k a = (k,k, —k, —k) (see figure [5]). The two Z 2 
symmetries of the quiver, one acting by interchanging nodes 1 -H- 4 and 2 -H- 3 and the other 
by interchanging nodes 1 -H- 2 and 3 -H- 4, allow us to set 



A 



2,j 



A 



,1 ■ 



A 



3,j 



A 



4,j 



(43) 



Moreover, in the saddle-point equations (34) it is consistent to set Xj = Xj as in the ABJM 



case, which reduces our task to finding a single eigenvalue distribution Aj. In passing to the 
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ki = k 



< > 



ko = k 



< > 



ki 



-k 



-k 



Figure 4: Four-node quiver diagram obtained as a particular case of the general quivers 
presented in figure |3j 

continuum limit, we should therefore set 



Vi = V2 



-2/4 = 2/ 



(44) 



The free energy functional (39) then becomes 
2kN 3 ' 2 



F[p,y] 



7T 



dxxp{x)y(x) + N 3/2 / dxp{x) 2 [n 2 + f{2y(x))] + o{N 3/2 ) . (45) 



In the paragraph following eq. (39) we discussed how for arbitrary p-node quivers we 



would not be able to solve for the y a themselves, but only for differences of consecutive y a , 
because the leading large N contribution to the free energy functional is invariant under the 
shifts y a — >■ y a + 5y for any function 5y. In the case of the (k, k, —k, —k) quiver we will, 
however, be able to determine the location of the eigenvalues exactly, because the ansatz 



(44) breaks this shift symmetry. 



In order to find the saddle points of (45) in the set (14), we should add a Lagrange 



multiplier /i to enforce the normalization condition for p and extremize the functional 

N 3 ' 2 ( f 
F[p, y} = F - { dx p(x) - 



(46) 



Let us assume the eigenvalue distribution is symmetric around x = y = and ranges 
between [— x*,x*\. Let us focus on the region where x > 0. Solving the equations of motion 
we obtain 



p(x) 



P 



y{x) 



2kn 2 x 
P 



if \y(x)\ < 



71 



(47) 
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Since p{x) > in this region, we have p > and y(x) > 0. Assuming y{x 3f ) < ir/2, we can 
find p in terms of x* from the normalization condition for p, and then express F in terms 
of and extremize it. The extremization yields x* = 2 1//4 7t/a/A; and y(x*) = Tv/y/2 > 7r/2, 
which suggests that the assumption y(x*) < n/2 might be wrong. One could imagine that 
y(x*) > 7r/2, but solving the saddle-point equations in the region where y > 7r/2 would yield 
p(x) < 0. 

The correct answer is y(x*) = ir/2, and in fact y(x) = n/2 on some interval [x 7T / 2 ,x 3f ] 
with < x n /2 < £*. On this interval, 



p{x) 



p — 2knx 

47T 3 



y( x ) = o > 



(48) 



where in obtaining these equations we only varied (46) with respect to p. The quantity x n /2 
can be obtained from setting j/(x 7r / 2 ) = vr/2 in (47): 



X-k/2 



JJ_ 

4nk 



(49) 



One can now find p by imposing the normalization condition for p, and then express the free 
energy F in terms of x* and extremize with respect to x*. The result is that 



x* = 2x^/2 = 2ir\! ^ , 




p = 47T' 



2k 




(50) 



The density of eigenvalues is constant on [— x n / 2 , 2^/2] an d then drops linearly to zero on 



Af=100, k={l,l -1 -1} 

y 







\r/4 




-4 -2 / 


2 4 







Af=100, k={\,\ -1 -1} 
P 

0.14 t 



I 0.10 

0.08 
0.06 
0.04 
0.02 



A 



Figure 5: Comparison between numerics and analytical prediction for the four- node quiver 
with k = {1, 1,-1,-1}. The dotted black lines represent the large N analytical prediction, 
and the orange dots represent numerical results. 



— x*, — 2^/2] and [x n /2, x*]. See figure p^for a comparison of this analytical prediction with a 
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numerical solution of the saddle-point equations. 

The free energy for this model can be computed from (45): 




7T 



^1/2^3/2 + o(jv 3/2 ) 



(51) 



Using (|40j), we infer that the gravity dual of the Chern-Simons quiver gauge theory with CS 
levels {k, k, —k, —k) is AdS^ x Y where the volume of the compact space Y is 



Vol(F) 



7T 

16k 



(52) 



Satisfyingly, this result is in agreement with the calculation of the corresponding integral 
representation given in [20] for k = 1, which we will review in section 4. 



k\ = k 



ko = k 



:: 



ki = —k k 3 = —k 

Figure 6: The Chern-Simons quiver gauge theory dual to AdS^ x Q 2 > 2 ' 2 /Z k as proposed 



m 



34,35 



Let us also note that this volume is the same as that of a orbifold of the Sasaki- 
Einstein space Q 2,2 ' 2 , which in turn is a Z2 orbifold of the coset space SU(2) x SU(2) x 
577(2)/ (U(l) x U(l)). If we denote the generators of the three SU(2) factors by J a, Jb, 
and Jc, then the two U(l) groups we are modding out by are generated by Ja3 + Jb3 and 
Ja3 + Jc3- Q 2 ' 2 ' 2 admits a toric Sasaki-Einstein metric, and a proposal for the Chern-Simons 
quiver gauge theory dual to AdS± x Q 2 ' 2 ' 2 /Zk was made in 34 35 . This proposal is quite 
similar to the (k, k, —k, —k) nonchiral quiver in figure El except it is chiral — see figure [6j 



Because of the chiral nature of the quiver, the corresponding matrix model that follows 
from [7] is somewhat different. Its analysis is beyond the scope of this paper. 
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3.4 Extremization of the Free Energy Functional and Symmetries 



Since the free energy functional (39) depends only on differences between consecutive y a , we 



find it convenient to introduce the notation Sy a = y a -\ — Ua an d to write k a = q a+ i — q a as 



in eq. (32). Equation (39) becomes 



n 3 / 2 r p n 3/2 r p 

F [P^Va] = -o— / dxxp(x)22q a 6y a {x) + — / dx p{xf ^ f(6y a (x)) + o(iV 3 / 2 ) . (53) 

a=l J a=l 

This expression should be extremized over the set 

C = ^(p, 8y a ) : J dx p{x) = 1; p(x) > and '^^5y a (x) = pointwisej' . (54) 

Since Y7 a =i ^Va = ®-> one cou ld either use this constraint to solve for one of the Sy a and 
extremize ( |53| ) only with respect to the remaining ones, or, as we will do, one could introduce 
a Lagrange multiplier u(x) that enforces the constraint and treat all 5y a on equal footing. 



Because of the normalization constraint (38) we also need a Lagrange multiplier p. We 
therefore will extremize 



N 3/2 ( f \ N 3 ' 2 , x ^ 

F\p,Sy a ] = F[p,5y a ] — p ( / dxp(x) - 1 J - J dx p(x)u(x) 2_^5y a (x) (55) 

~ " _ a=l 



2tt 



instead of (53). Suppose a saddle point exists. As long as p(x) > 0, the saddle-point 



eigenvalue distribution should satisfy the equations 



^2 l 27T f( 5 ya( X ))p( X ) + (laX - V{x)) 5y a {x)\ = (I , 

%f(5y a (x))p(x) + q a x = u(x) 



a=l 



(56a) 
(56b) 



The extremization problem has the following discrete symmetries: 



The free energy functional (53) has a Z 2 symmetry under which q a and 5y a all change 



sign, so in the large N limit the partition function and the free energy are also invariant 
under sending q a — > —q a for all a. This symmetry acts as k a — > —k a and is therefore a 
parity transformation. 



Equation (53) is invariant under an overall shift of all the q a . This symmetry was to 



be expected given that, after all, the original integral (33) depends only on k a , which 
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are differences of consecutive q a . 

Interestingly, the free energy functional we are extremizing is invariant under per- 
mutations of the q a and 8y a , so the partition function and the free energy will also 
be invariant under permutations of the q a . Up to order 0(N°) shifts in the ranks 
of the gauge groups, which should be dropped in the large N limit we are taking, 
such permutations correspond to Seiberg dualities in the M = 2 Chern-Simons gauge 



theories 21-23 . 



Some of the symmetries discussed above correspond to the action of the dihedral group 
D p on the CS levels k a . Our formalism shows that to leading order in N the free energy is 
in fact invariant under a larger symmetry group that acts on the q a and that includes the 
dihedral group. 



3.5 Three-Node Quivers 

Let us now compute the free energy for arbitrary three- node quivers with CS levels (ki,k 2 ,k 3 ) 
satisfying k\ + k 2 + k 3 = 0. Since the k a sum to zero, two of them must have the same sign 
and be smaller in absolute value than the third. Let us begin by studying the particular case 
where k 2 > ki > and k 3 < 0. For simplicity, we choose Yl!=i = 0, which implies 



qi 



2fci + k 2 



Q2 



h - k 2 



93 



h + 2k 2 



(57) 



and we have g 3 > > q 2 > q\ and |g 3 | > > \q 2 \. The solution to eqs. (56) is symmetric 
about x = 5y a = 0, and when x > it breaks into three regions: 



< x < 



37r<? 3 



< X < 



37r 2 xq a 



6tt 3 ' 



3%q 3 

» . 

37rgi ' 

y 4vrp 2 ' y 



(92 



Airp 



qi)X _ 7T 

2 



P 



2/i — 37rg 3 a; 

67T 3 



(58a) 



(58b) 



6y 3 = it 



(58c) 
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37rgi 



< x < 



K(q 3 ~ qi) ' 

5yi = -7T , 6y 2 = , 5y 3 = vr . 
M + (ft - 93)^ 



(58d) 



2tt 3 



The first region ends when one of the three differences Sy a reaches ±7r. The relations 



Af=100, k={l,3, -4} 




jV=100, ^={1,3,-4} 
P 



0.15 
0.10 
0.05 



1 



Figure 7: Comparison between numerics and analytical results for a three-node quiver. The 
dotted black lines represent the analytical large N approximation, while the orange dots 
represent numerical results. 

between the q a imply that at the end of the first region Sy 3 = n, while \Syi\ = k \ qi \ / \ q 3 \ < tt 
and \5y 2 \ — 7r \q 2 \ / | ^3 1 < 7r. Throughout the second region 8y 3 = ir. The second region ends 
when Syi or Sy 2 reaches ±n. When q\ = q 2 , the third region is absent. When qi < q 2 < 0, 
in this region 5y 2 is monotonically increasing and Syi is monotonically decreasing, and since 
^2^Ua — it must be that Syi reaches —it next. In the third region the 5y a are all constant 
and the density p decreases linearly to zero. See figure [7] for a particular example. 
The normalization condition on p yields 



P = 7T 



18gi(gi + g 2 )(2gi + q 2 ) 
q\ ~ 5g? - 5gig 2 



l 2{k x + k 2 ){k 2 -k 3 ){k x -k 3 ) 
(kik 2 - k x k 3 - k 2 k 3 ) 



(59) 



Performing the integral (39), one obtains 



N 3 / 2 p N 3 / 2 nV2 {h + hXh-hXh-h) 



3tt 



k\k 2 - k x k 3 - k 2 k 3 



(60) 



Given the free energy in the case k 3 < < ki < k 2 , it is actually possible to compute 
the free energy for any three-node quivers. Indeed, since in the case where there are only 
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three nodes a permutation of the k a can be thought of as a relabeling of the nodes, the free 
energy must be invariant under all such permutations. In addition, the free energy must be 
invariant under sending k a — > —k a according to the second discrete symmetry discussed at 



the end of section |3.4[ Combining these two properties, one can find the free energy of an 
arbitrary quiver with CS levels k a by constructing the new CS levels k\ = mindfexl , \k 2 \ , \ k s \), 
ks = — max(|&i| , l&^l , |*3|)> and k 2 = ~ *i — k% that satisfy k% < < k\ < k2 and for which 



eq. (60) holds. The unique extension of (60) that gives the correct answer for arbitrary CS 



levels is 



(1*1 


+ 


|* 2 


)(N 


+ 


\k 


s|)( 


h 


+ 


h 


) 




fell 


1*2 


+ 


*i 


| | 


k 3 


+ 


k 2 \ 


1*3 





(61) 



3 

Quite remarkably, this formula, whose derivation is based solely on gauge theory argu- 



ments, agrees with the supergravity prediction: Using (40), one can reproduce the volume 
of a Z gcd { fcl fc2 j C3 } orbifold of a compact Eschenburg space. The Eschenburg space is specified 
by three relatively prime integers t a , and its volume is fl9] 

(62) 



Vo\{S{t u t 2i t 3 )) hh + Uh + hh 



Vol(,S 7 ) (tl+* 2 )(<2+*3)(*l+*3)' 

In terms of the k a , the positive integers t a are t a — \k a \ / gcd{fci, k 2 ,k 3 } [13], so 



Vol(F) 1 Vol^Ma,^)) 


*i| 


\k 2 


+ \> 


dl 


1*: 


.1 + 


k 2 \ 


1*3 




Vol(^) gcd^,^} Vol(^) (\h 


+ 


\k 2 


m 


1 + 1' 


*|)( 


*i| 


+ N) 



(63) 



in agreement with (40) and (61). 



3.6 General Four-Node Quivers 

We can also compute the leading large N contribution to the free energy for arbitrary four- 
node quivers. Let us first examine the case where g 4 > q 2 > q\ > q 3 and \q±\ is the largest 
among the q a . It is convenient to require Ylt=i la = since many of the intermediate formulae 
simplify under this assumption. Then we have q^ > > q\ > q% and {q^l > \q%\ > \q\\ > \q 2 \. 



As in the three-node case, the solution to eqs. (56) is symmetric about x = Sy a = 0, and 
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when x > it breaks into three regions: 



< x < 



4vrg 4 



< x < 



47rg 4 
47rg 3 



8tt 3 



(64a) 



(64b) 



%3 



4vrg 3 



< a; < 



27r(g 2 + 94) 

Syi 



(3gi + g 4 )x _ 7T 

67rp 3 ' 

(3g 3 + q 4 )x _ 7T 
67rp 3 ' 

(gl - ^2)2; 



^2/2 



(3g 2 + ^4)2; 7r 
67rp 3 



5y4 = ir . 



P 



3/i — inq^x 
167^ 



(64c) 



hi 



$V3 



Sy4 = 7T 



(gg - qi)x 

_ + (gg ~ g4)7Tg 

P " 4vr 3 



The first region ends where £7/4 reaches 7r. At this endpoint \Sy a \ = ir \q a \ / \q±\ < it for 



N=45, £={3,-4,7,-6} 
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Figure 8: Comparison between numerics and analytical results for a four-node quiver. The 
dotted black lines represent the analytical large N approximation, while the orange dots 
represent numerical results. 

a = 1,2,3. The second region ends where Sy 3 = —ir. At this endpoint 5y\ = n(qi — 
g2)/(gi + I2 — 2g 3 ), and since qi > g 3 and g 2 > g 3 , by the triangle inequality it follows that 
\qi — 92 1 < qi + g2 — 2g 3 , so \Sy\\ < n. Similarly, \5y 2 \ < n also. Lastly, if g 2 = gi, the third 
region does not exist. When g 2 > gi and gi < 0, is monotonically decreasing and 5y 2 
is monotonically increasing in the third region, and this region ends where 5yi = —71 and 
5y 2 = 7r. See figure [8] for an example. 
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From J dx p(x) = 1, one can find that /i is given by 



8^ = /I_I + 4ta + *) + ^2__ (65) 

A* V 93 94 (92 + 94) 2 92 + 94 



The free energy is 



iV 3 / 2 /i 8vriV 3 / 2 / 1 1 4(g 2 + 9s) 12 



+ W2 - " + • (66) 



3tt 3 \q 3 g 4 (g 2 + 94) 2 92 + 94 



-1/2 



to 



Given eq. (66), one can use the symmetries we discussed at the end of section 3.4 
compute the free energy of a quiver gauge theory with arbitrary q a . Indeed, one can define 
q a to be a permutation of the four numbers q a — \ Ylt=\ 9fc that gi ves 1 94 1 > 1 93 1 > |<7i| > 1 92 1- 
If g 4 is negative, one should flip the sign of all q a , so we can assume g 4 > 0. By construction, 
the q a sum to zero, so the second and third largest in absolute value, namely q 3 and q~i, are 



negative. Therefore, the q a satisfy all the assumptions under which eq. (66) was derived, and 
since the free energy does not change in going from q a to q a , one can plug the q a into eq. (66) 
to find the free energy of an arbitrary four-node quiver theory. The unique extension of (66) 
to arbitrary q a can also be written as 



p _ N 3/2 nV2 nti (Etx M j (67) 

3 V ^(a,b)^(c,d)^(e,f) \lab\ \Qcd\ \lef \ ~ ZJ(o,6,c) \1ab\ \Qbc\ \Qca\ 

where q a j, denotes q a — qb-, and in the denominator the first sum is over distinct unordered 



pairs of numbers from 1 to 4 while the second sum is over unordered triplets. Using eq. (40) 
we obtain a prediction for the volume of the compact space Y: 



Vbl(y) _ S(o,6)^(c,d)^(e,/) \lab\ \lcd\ \q e f \ ~ Y,(a,b,c) l ?a6 l l ?bc l l ?ca l 

VoiM" n:=i (Eti M) 

4 A General Formula and its Tests 



(68) 



Equations (67) and (68) suggest a generalization to arbitrary p-node quivers. Note first that 



the numerator of eq. (68) is a sum over all possible graphs with 4 nodes and 3 edges from 
which we subtract the sum over all cyclic graphs with 4 nodes and 3 edges, yielding a sum 
over all possible trees. 

We conjecture that for a p-node quiver, the volume of the tri-Sasaki Einstein space Y 



24 



(normalized so that R mn = 6g mn ) is given by 

Vol(y) _ Y,(V,E)eTH(a,b)eE 1 9a - 76 1 



Vol(S 7 ) 



ELl 1 7a - 7b I 



(69) 



where T is the set of all trees (acyclic connected graphs) with nodes V = {1, 2, . . . ,p} and 

edges 



£ = {(ai, fei), (a 2 , 6 2 ), • • • , (ap-i, V-0) • 



(70) 



A standard result in graph theory states that trees with p nodes have p — 1 edges. 



The conjecture in eq. (69) is consistent with the results from two-, three-, and four-node 



quivers, and we also checked it for five- and six-node quivers. This formula is invariant 



under all the symmetries discussed at the end of section 3.4 In particular, a quite nontrivial 
check of our approach is that this formula is invariant under the Seiberg dualities described 
in 21 23 . The connection we observe between large N matrix integrals and sums over the 



tree graphs is reminiscent of the connection between matrix models for 2-d quantum gravity 



and the Kontsevich matrix model which generates ribbon graphs 36 



An integral representation of volumes of tri-Sasaki Einstein spaces was given by Yee [20]. 



In general, our spaces Y are orbifolds of those considered in 20 , where k = gcd{k a }. To 



simplify the following discussion, let us focus on the k — 1 case. In this case [20], 



Vol(Y) 



2P~ 2 7T 4 



3Vol([/(l) 



(71) 



Here, Vol (U(l) p ~ 2 ) is the volume of a unit cell in the (p — 2)-dimensional lattice defined by 
the identifications ~ + rjj, where rjj satisfy Y^Zi QiVj 2vrZ for all a = 1, . . . ,p. The 
Q J a span the kernel of 

( i i 1 ... i \ 

(72) 

(The Q J a are taken to be relatively prime here.) In the U(l) p Chern-Simons gauge theory, the 
Q J a are the charges of the bifundamental fields under the unbroken U(l) p ~ 2 symmetry [j~3~] . 
We can take a spanning set of Qi to be, for a fixed j, Q{ = q<i — qj, Q 3 2 = qj — qi, and 
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Qj = qi — q2 with all other Q{ = 0. For this choice of Q 3 a , the volume of U{l) p 2 is 



Vol (f/(ir 2 ) = (27r)P " 3 . (73) 
\Qi - gal 

Note that Vol(Y) is invariant under permutation of the g . Although we have not carried 
out the integral in general, we can investigate specific cases with ease. For example, for the 



choice q = (3, 2, 1, 2), corresponding to the k = (1, 1, —1, —1) quiver, both our formulae (69) 



and (71) give Vol(Y) = 7T /16. A more nontrivial choice is q = (3,2,1,5) for which both 
formulae yield 1397i~ 4 /4725. By evaluating (71) numerically, we were able to check agreement 



with (69) in a number of randomly selected cases for p = 4, 5, and 6. 



The volume formula (69) is invariant under a shift q a — > q a + 1. In the type IIB brane 
construction, which involves a sequence of (1, q a ) 5-branes, this symmetry corresponds to the 
T transformation of the SX(2,Z) S-duality group. We could use the SL(2, Z) symmetry to 
generalize the free energy to theories whose brane constructions involve general u a = (p a , qa) 
5-branes. This generalization is accomplished by replacing the differences \q a — qb\ in the 
volume formula with \u a A Ub\ = \p a qb — Pbq a \- For special cases where some of the p a vanish, 
this formula describes theories with fields in the fundamental representation. For example, 
for the ABJM model with Nf flavors, corresponding to u\ = (1, k),U2 = (1, 0), u% = (0, Nf), 
our formula predicts 

vag) = 2k + N f (7A) 

Vol(S 7 ) 2(k + N f y 1 1 



This equation agrees with the explicit matrix model calculation 37 and with the volumes 
of Eschenburg spaces S(Nf, Nf, k) fl~9]. 



5 Discussion 

In this paper we have studied p-matrix models describing certain U(N) P Chern-Simons quiver 
gauge theories with M = 3 supersymmetry. In the large N limit these theories are dual 
to eleven-dimensional supergravity on AdS^ x Y, where Y is a tri-Sasaki Einstein space. 
By finding an analytical large iV limit of the matrix integrals, we were able to check the 
supergravity prediction that the logarithm of the partition function of the gauge theories 
on S 3 should grow as N 3 ^ 2 . In AdS^ x Y the coefficient of proportionality depends on the 
volume of the compact spaces Y, so we could compare our gauge theory results with the 
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volumes computed earlier using geometric techniques 19,20 . These successful comparisons 



constitute new detailed tests of the AdS 4 /CFT 3 dualities. In eq. (69) we conjectured an 
explicit combinatorial volume formula for arbitrary p. It should be possible to derive this 



formula in an independent way using algebraic geometry techniques similar to those in 38 



Quite generally, the main difficulty in solving matrix models is that the interactions 
between the eigenvalues are long-ranged, and the saddle-point approximation yields integral 
equations in the continuum limit. Remarkably, in solving the models described in this paper, 
one can set up an approximation scheme where the eigenvalue distributions can be found 
by solving algebraic equations. The limit in which the saddle-point equations simplify is 
the limit of "large cuts" where the eigenvalues grow as an appropriate positive power of 
N. Perhaps the key insight in solving these matrix models was that the long-range forces 
between the eigenvalues can be made to vanish by choosing the distribution of the real parts 
of the eigenvalues to be the same for each set of eigenvalues. The remaining interaction forces 
between the eigenvalues are short-ranged, and that is the reason why in the right variables 
the saddle-point equations were local and algebraic in the large N limit. 

While we worked in the limit where N is sent to infinity and the Chern-Simons levels 
k a are kept fixed, it is of obvious further interest to relax these assumptions and study 1/N 
corrections. In doing so, a subtle issue that needs a better understanding is the imaginary 
part of the free energy. At first sight, the imaginary part in the ABJM model is of order 
O(N). On the other hand, one could argue that this imaginary part is only defined modulo 
27r because a shift of the free energy by an integer multiple of 27rz leaves the partition function 
unchanged. 

Another interesting generalization of our results is to solve the matrix model in the scaling 
limit where the Chern-Simons levels are sent to infinity, with N/k a kept finite. One could 
calculate the free energy as a function of the 't Hooft-like couplings N/k a and check that, as 
predicted by the AdS/CFT correspondence, it should interpolate between an N 2 behavior 
at small N/k a dictated by perturbation theory and the fc 1//2 iV 3//2 behavior at large N/k a that 
we found. For p = 2 this check was performed in [5] by computing the resolvent of the 
matrix model using the techniques developed in 10 . We believe a similar check should also 
be possible for the Af = 3 theories studied in this paper, using perhaps similar techniques. 
Such an approach should also provide access to the ABJ-like cases where the ranks of the p 
gauge groups are not equal. 

Finally, it would be interesting to investigate whether the large N matrix integrals we 
have calculated play a role in four-dimensional gauge theories, for example, in the 4-d "parent 



27 



theories" 
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of the 3-d Chern-Simons models we have studied. 
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A Taking the Continuum Limit of the Free Energy 



We start with the matrix integral eq. (|33|) and write it as 

- f(x ^ n ^ , (75) 



e 



where we divide the free energy into the following three pieces: F = F ext + F int + -F const We 
have defined F ext to be the contribution to the free energy from the external potential 

a,i 

The contribution to F from the eigenvalue interactions is 



v H 4sinh 2 (^=^Y 

io g n— — h V (77) 

a A in M 2cosh(^±M A 



3 In the following, we will use the product formula \og(xy) = \ogx + logy, although this is strictly 
speaking only correct up to integer multiples of 2ni. The extra contributions would not affect the saddle- 
point equations. 
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Finally, there is an overall normalization F const chosen to be consistent with |5l: 



i^const =P (log N\ + N log 2vr) 



(78 



It will turn out that -F CO nst does not contribute at leading order in our large N expansion. 

To take the continuum limit, we begin with the assumptions justified in the text of the 
paper that the eigenvalue distributions lie along curves symmetric about the origin in the 
complex A plane, that we can order the eigenvalues such that Re A a j = Re A^j for any a and 



b, and that |ReA aj -| 3> |Im A a j|. More specifically we assume the following variant of (35) 



K,j = N a Xj + iy a7 



3 i 



(79) 



where a > 0. 

Taking the continuum limit of F ext is straightforward. The leading term in iV cancels 
because k a = 0, and we are left with 



N a 
2tT 



Fad = ~— ^2 koPajVaj + °(N) 



(80) 



a- J 



Letting the real parts of the eigenvalue distributions extend from — x* to x*, we can introduce 
an eigenvalue density p(x) and approximate the sum over j as an integral over x: 

N 1+a k f x * 

F ext = J2 k * / xy a (x)p(x) dx + 0(N). (81) 

a 



Taking the continuum limit of F int is more involved. We begin by reorganizing the 
products: 



int 




_|_ e ~K,i+>>a+i,j}(l _j_ e -Aa,i+A _i,j) J pj 2 cosh ^ Xa ' i ~ Xa + 1 ' 



(82) 
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We then convert the logarithm of the product into a sum over logarithms: 



fl- = E 

a=l 



i>j n=l 



+ 



+ log I 2 cosh 



A a ,i — -^a+l,i 



(83) 



Making use of the assumption (79) and taking the continuum limit, the interaction energy 
reduces to 



fl- = E 



a=l 



iV log 2 cos 



' 00 _/\r 2 



+ / ^ — 2e ( - Aa(:c)+Aa(2: ' )) "+ 



n=l 



-( — l) n ^ e (- A a(a:)+Aa+i(a;'))n _|_ g(-A a (a;)+A a _i(a:'))n j j p( x ') 



(84) 



p(x)dx . 



We now estimate the integral over x' in the above expression for F- mt - Consider the 
following related integral: 



1 



^ Xb ^ +X ^ n p(x')dx' 



n 
1 dx 
n dX n 



i-X h (x)+X a (x'))n dx ' / a 



p(x)e 



(-\ b {x)+\ a (x))n 



1 r 



,(-A t ( a Q+A,(a'))n d „f„/' 



p(x') dx' (85) 



+ 



Given (79), the integral in the second line will be suppressed by a factor of 1/N a compared 



with the boundary term. The boundary contribution from x* will be suppressed by an 
exponential amount because Re [A (— x*)] < Re[A&(x)]. The last line of (85) reduces to 



71/ 

Introducing the notation 8y a = y a -\ — y a , the interaction energy reduces to 

„ x* 

p(xf 



a=l n=l 

—x, 



F- mt = N 2 ~ a J2 / Yl 2 ~ ( _1 ) n {e' in5ya+l{x) + e inSya{x) ) 



n- 



dx + 0(N 2 ~ 2a ,N). (87) 



We are tacitly assuming that a < 1 and so can drop the order N term from the energy. 
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Reorganizing the sum over a, we can write this energy as 



int 



a=l 



£ / f(6y a )p(x) 2 dx + 0(N 2 - 2a ,N). 



where we have defined the function 



/(y)=£^[l-(-ircosn2,] 



39) 



n=l 



Clearly / is a periodic function of y with period 27r. Recall the Fourier series expansion for 
y 2 in the domain — tt < y < tt: 



00 4f_iy 
2/ 2 = £^^cosm/ + 2C(2) 



n=l 



(90) 



In the fundamental domain — tt < y < it, the function / is thus 



(91) 



B A More Detailed Check for ABJM Theory 



We explore solutions to eqs. (17) where \y\ > tt/2. Based on the numerical results in 



section 2.2 we expect the eigenvalue distributions to be invariant under Aj — > — A.; and 



Aj —> — Aj, which implies y(—x) = —y(x) and p(—x) = p(x). Assuming k > 0, we will focus 



only on the x > region. Plugging (12) into (|17|) one obtains 



7T 2 kx 



2{i 



if |y(x)|< 



TT 



(92) 



and 



/c(/x — 2kx7r) 



y(x) 



7r(2/i — 3/cx7r) 



„ 7T 37T 

if 2 - y( x > - y 



(93) 



47T 3 ' ' 2p — Akxn 

and so on. From eq. (|9~2| we infer that /x > and y(x) > if x > 0. We could have in 



principle also allowed y(x) = n/2 over some range of x, but then the first equation in (17) 
would imply that x = fi/nk, so y(x) could equal tt/2 only on a set of measure zero. 

Assuming a connected distribution of eigenvalues of each type where p is supported on 
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[— x*,x*] for some > 0, there are two possibilities: Either y{x*) > n/2 or y(x*) < tt/2. 
Assuming y(x*) > 7r/2 we immediately reach a contradiction. Indeed, consider the point 
£71-/2 = fi/itk where y(x n / 2 ) = vr/2 and eq. (92) joins onto eq. (93). For x > x^/2, eq. (93) 



implies p(x) < 0, which contradicts the assumption that p(x) > 0. It must be that p(x) 
for x > x n / 2 and thus \y(x)\ < n/2 for our eigenvalue distribution. 
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